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ABSTRACT 


Aims. We investigate the electron acceleration in convective electric fields of cascading magnetic 
reconnection in a fiaring solar corona and show the resulting hard X-ray (HXR) radiation spectra 
caused by Bremsstrahlung for the coronal source. 

Methods. We perform test particle calculation of electron motions in the framework of a guiding 
center approximation. The electromagnetic fields and their derivatives along electron trajectories 
are obtained by linearly interpolating the results of high-resolution adaptive mesh refinement 
(AMR) MHD simulations of cascading magnetic reconnection. Hard X-ray (HXR) spectra are 
calculated using an optically thin Bremsstrahlung model. 

Results. Magnetic gradients and curvatures in cascading reconnection current sheet accelerate 
electrons: trapped in magnetic islands, precipitating to the chromosphere and ejected into the 
interplanetary space. The final location of an electron is determined by its initial position, pitch 
angle and velocity. These initial conditions also influence electron acceleration efficiency. Most 
of electrons have enhanced perpendicular energy. Magnetic curvature and gradient driven ac¬ 
celeration efficiency along the magnetic field increases with the increase of the magnetic field 
resolution. It was shown that for a sufficiently high resolution of the electromagnetic fields, adap¬ 
tive mesh refinement (AMR) MHD simulations are needed. Trapped electrons can be accelerated 
to energies up to half of MeV, precipitating ones to more than 60 keV in the refined smaller-scale 
magnetic structures. Trapped electrons are considered to cause the observed bright spots along 
coronal mass ejection CME-trailing current sheets as well as the flare loop-top HXR emissions. 
Another observable effect is the locations of the precipitating energetic electrons with respect 
to the polarity inversion line (PIL). A stronger asymmetry of the energetic electron precipita- 
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tions around the PIL is obtained with sufficiently well resolved smaller-scale magnetic structures 
than that in the coarsely resolved ones due to the asymmetry in the parallel magnetic curvature 
acceleration. 

Key words. Solar flare electron acceleration - Adiabatic acceleration - Hard X-ray - UV radiation 
- Cascading reconnection 


1. Introduction 


Since the first recorded white light observations of solar flares ( |Carrington| 1 859t |Hodgson| 1 859| ), 
sophisticated ground-based and space-born solar techniques have been introduced to investigate 
the physics of the sun. Recently space telescopes like SoHO, Yokoh, RHESSI, Hinode and SDO, 
have revealed many detailed observations covering broad wavelength ranges at a high temporal, 
spatial and spectral resolution. 

Generally, it is accepted that the energy of solar flares comes from stressed, non-potential, 
current-carrying coronal magnetic fields being released by magnetic reconnection. About 10 to 


50% of the flare energy may be transferred to energetic electrons and ions (e.g. Lin & Hudson 


1976| ). In some cases energetic electrons alone carried away 50% of the flare energy (e.g. Miller 


et al.|1997) ), being accelerated to energies up to 10 - 100 MeV (e.g. | Aschwanden|2002| ) . 


The prime diagnostic of accelerated electrons in solar flares is the HXR radiation they cause. 
Two main components were identified in HXR light curves : a sharply increasing component and a 
slowly varying one. The sharp increase happens within 0.5-5 s after the initial flaring (e.g. |Holman| 
et al.||20lT[ [Zharkova et al.||20TT] ). This indicates that within sub-second electrons are locally 


accelerated in excess of a few MeV. The slowly varying component lasts as long as flares continue, 
i.e., electron energization continues. 

Using high-resolution imaging, HXR location in the sun has been solved with a few arcseconds 
resolution. Solar observations have shown HXR emissions from the foot-points of flaring coronal 
structures. Recently based on the classical CSHKP (see [Priest & Forbes|2002| for a review) solar 
flare reconnection model and solar flare observations near the limb of the solar disk, HXRs were 
found also at the flare loop tops (e.g., [Masuda et al.|1994t[Gordovskyy et al.|2010b[ with Yohkoh 
observation). 

Although a substantial progress was made in observations, it is still an open question by which 
mechanisms the flare electrons are accelerated. Mostly suggested mechanisms can be divided 
into three classes (1) acceleration by direct current (DC) electric fields (see, e.g., [Zharkova~& 


[Gordovsky^[2004| [2005a|b[ ), (2) stochastic acceleration (see, e.g., [Vlahos & Cargill[[2009[ ) and 
(3) shock acceleration (see, e.g., [Aschwanden[[2002t [Benz[[200^ . Observations show also that 
different flares produce different HXR spectra changing with time and their locations with respect 
to the polarity inversion line (PIL) (e.g., [Zharkova et JI][2011| ). All these features can hardly be 
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explained by one single acceleration mechanism. Therefore flare energetic particles are perhaps 
accelerated by different mechanisms at different time and in different places while the flares last. 

In order to validate acceleration mechanism, it is appropriate to carry out test particle calcula¬ 
tions. The electron acceleration in the vicinity of a single reconnection X-point, e.g., was investi¬ 


gated by [Zharkova & Gordovskyy|2004| |2005a|bt [Wood & Neukirch|2005[ [Priest & Titov|1996| ). 
Analytic reconnection models were used as well as the results of ideal or resistive MHD numerical 
simulations. [Martens & Young|1990| e.g., used MHD simulation results to study the particle mo¬ 
tions and accelerations in current sheets. In all these models the magnetic field parallel component 
of the DC electric field E = -u x B -h 77 J causes strong particle acceleration if only rj was chosen ap¬ 
propriately. However, the prescription of r] in the resistive MHD simulations is usually ad hoc and 
arbitrary. Meanwhile in the collisionless corona, the concept of collisional resistivity 77 is largely 
unapplicable. Microphysical effects have to be taken into account. [Silin et ar][2005[ [Buchner ^ 
Elkina[[2006[ e.g., have shown that considering possible micro-turbulence strong parallel electric 
fields must be confined in narrow channels of the ion inertia scale size (see also J. Buchner and 
W. Daughton 2007, section 3.5 in [Bim & Priest[ ( |2007[ )). Macroscopic MHD simulations, on the 
other hand, is better to be used to investigate the electron acceleration in the convective electric 
fields (E = -u X B). [Vekstein & Browning|1997[ and [Guo et al.|2010| e.g., analysed the particle 
acceleration in the convective electric fields around and at a magnetic null point, respectively. 


Vekstein & Browning] 1997 [ used an analytically prescribed magnetic field with an added uni¬ 


form electric field in the perpendicular direction to calculate the test particle guiding center motions 
near a reconnection X-point in a 2D geometry. They restricted the test particle orbits far away from 
the X-point since the guiding center approximation breaks down in a null-point. They considered 
particle parallel acceleration due to the E x B drift effects and neglected the effects of magnetic 
gradients and curvatures by launching only particles with very small initial parallel velocities and 
magnetic moments. As other authors before (e.g., [Burkhart et al.|1990| ), they found the final kinetic 
energy of the most accelerated particle is proportional to where E = -u x B. They also 
assessed the spectral index of the accelerated particles as being about 1.7, the corresponding HXR 
spectral index would be around 2.7 utilizing a simple relation = 6 + \ (where 6 is the electron 
spectral index and js stands for the index of emitted HXR spectrum) which is valid within the thin 
target model ( [Datlowe & Lin|197'3] ). 

Contrary to [Vekstein & Browning] 1997][Guo et al.|2010[ took the output of a 3D MHD simula¬ 
tion of magnetic null point reconnection to study the electron and proton acceleration at a 3D null 
point in the convective electric field. Every test particle is traced by solving the full equations of 
motion. This is necessary since the guiding center approximation breaks down at a magnetic null 
point. They investigated the influence of the convective speed on particle acceleration by rescaling 
it. They found that all particles are more efficiently accelerated with a larger convective speed. 
Particle energy can be up to energies of the order of 2 MeV (proton) and 3 keV (electron) from 
initial thermal energy of about 200 eV. The reason is that non-adiabatic (demagnetized) particles 
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can easily be accelerated in the direction perpendicular to the magnetic field. Particle have to un¬ 
dergo strong perpendicular drift to be substantially accelerated. But particle final parallel energy 
can still dominate its final total kinetic energy. Because of the much smaller gyroradius, electrons 
are demagnetized in smaller regions than the protons, protons are accelerated to higher energies 
than electrons. These authors also studied the infiuence of the initial energy on particle accelera¬ 
tion: higher initial energies lead protons to be stronger accelerated, while the final kinetic energy 
of electrons were not infiuenced essentially. 


In the studies of |Vekstein & Browning] 1997] and |Guo et al.|2010| there were only one magnetic 
X- or null point. [Krucker et al.||2008] claimed that electron DC-acceleration at only one recon¬ 
nection X-point (e.g., [Zharkova & Gordovskyy|[^Q5at [Wood & Neukirch]|2005| ) cannot explain 
the huge number of accelerated electrons inferred from HXR observations. The possibility of par¬ 
ticle acceleration by cascading reconnection was mentioned firstly by jShibata & Tanuma]|2QQl| 
They conjectured that magnetic islands could be formed at many scales by tearing-mode instabil¬ 
ities of the stretching current sheet. Later this concept are confirmed by theoretical approaches 
(e.g., |Loureiro et al.|2007[|Uzdensky et al.|2010| ), observations (e.g., |Hoshino et al.| 1994 t[Karlicky| 
2004 j ), AMR MHD simulations (e.g., |Barta et al.|20TT] ) and particle in cell (PIC) simulations (e.g.. 


Karlicky et al.|2012| . The electron acceleration by many reconnection sites was studied by [Li & 


|Lin|2012| and |Gordovskyy et al.|2010a|b| In their studies, however, they assumed arbitrary ad-hoc 
prescribed anomalous resistivity models to reveal the accelerating fields. As well as the number of 
X-points in their studies were obtained by periodically repeating the simulation domain. 


Only the particle acceleration in the convective electric fields E = -u x B, however, is inde¬ 
pendent on any ad hoc assumption about anomalous resistivity. In order to understand its possible 
acceleration effects, we use the results of AMR-MHD simulations of multiple island formations 
by cascading reconnection ( jBarta et al.|2010|[2011| ). Those simulations have shown that cascad¬ 
ing reconnection forms differently sized magnetic islands where electrons can be accelerated (see 
Sect|^. We use two different magnetic structure resolutions to investigated the resolution infiuence 
on electron accelerations. We studied the electron acceleration by cascading reconnection not only 
near the X-points but also in the magnetic islands in the framework of a guiding center approxima¬ 
tion ( [Northrop 1 1963 see Sectj^. The resulting HXR emissions by energetic electron non-thermal 
Bremsstrahlung are derived using a optically thin Bremsstrahlung method ( [Brown[1971t[Tandberg-[ 
Hanssen & Emslie|1988[ ) to compare with flare HXR observations. In Sect|^ electron acceleration 


dependence on initial conditions, different acceleration factors in the parallel direction, acceler¬ 
ation in different (parallel and perpendicular) direction, as well as trajectories, are investigated 
for trapped (Sect [4.1[ ) and precipitating (Sect [4.2[ ) electrons. Finally the results are discussed and 
conclusion are drawn in SectlS] 
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2. Electromagnetic fields of cascading reconnection 


In this study we aim to investigate the particle acceleration in the convective electric fields of 
differently resolved cascading reconnection current sheets trailing a flaring arcade behind an ejected 
flux rope (cf. |Lin & Forbes |2000| ). 

The fields of cascading magnetic reconnection are obtained by means of a 2.5D AMR MHD 
simulation ( [Barta et al.|2010| 2011| ). In traditional MHD simulations, there are only uniform grid 
points. Unfortunately the sub-grid physics become important when the current sheet width and the 
non-idea plasma domain become thinner than the numerical grid size. Hence, traditional coarse 
MHD simulation cannot study smaller-scale processes of anticipated cascading reconnection. In 
order to resolve smaller-scale magnetic structures in the thinner current sheets, we use simulation 
results which can cover an as large as possible scale range. The high resolution AMR MHD 
technique allows the description of smaller-scale magnetic structures. For that sake the refined 
mesh is used when the current sheet width becomes comparable with the initial coarse grid size. 

The AMR algorithm works as follows: If at the time-step t + lS.t some coarse grids are detected 
containing thin current sheet, then they will locally be split into sub-boxes with 10x10 grid-points 
in the sub-system. After such refined meshes are initialized, the necessary more detailed plasma 
and field values are obtained by interpolating their parent coarse system values at the last time step 
(t). Then the dynamics of both the newly created and the pre-existing refined meshes are evolved 
in time {t ^ t + At) with an accordingly refined time-step. After that the plasma and field values 
at the parent coarse mesh are replaced by averaging the quantities obtained from its corresponding 
refined meshes at time-step t At. The influence of the global dynamics on the refined meshes 
are considered by interpolating boundary conditions in time and space. This refinement is repeated 
until the whole simulation is over (see [Barta et al.|2010| . 

So there are two sets of electromagnetic field data obtained by the AMR MHD simulation: 
one for a simulation on the coarse meshes alone and another with the refined meshes which pro¬ 
vides even smaller-scale structures of magnetic fields (see Fig|^. The MHD simulation results 
are restricted to 2.5D, i.e. two dimensional geometry but three dimensional plasma velocities and 
magnetic fields. This assumption is reasonable since observations have shown that the extended 
solar flare arcades typically having much larger extend along the polarity-inversion line (PIT) than 
across the PIT. 

The coordinate system is shown in Figj^ the x and y-axis are directed along and perpendicular 
to the current sheet, respectively. The current sheet center is located at y = 0, while the z-axis 
is pointing along the PIT located at (v = 0, y = 0). In this direction, every value is invariant 
i.e. didz = 0. The coarse resolution contains 6400 x 800 points in the vertical (x-axis) and right 
half of horizontal (positive y-axis) direction. A mirroring boundary is used at y = 0 for the left 
half box: p, Ux, By, and U are symmetric while Uy, B^ are anti-symmetric. For the upper 
and right sides, free boundary conditions are used: all quantities should satisfy the von Neumann 
prescription d/dn = 0 except of the normal magnetic field B^ and the total energy density U. Bn 
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and U are used to fulfill V • B = 0. At the bottom, a symmetric boundary condition {Q{-y) = Q{y)) 
is used for p, U and the anti-symmetric relation Q(-y) = -Q(y) is assumed for By. The 

plasma is always static u = 0 at the bottom. 


A generalized Harris-type current sheet is chosen as the initial state of the AMR MHD simula¬ 
tion ( |Barta et al.|2010[[2011| ): 


A(x,y,z',t = 0) = - 5 x 0 In 
Bz(x,y,zU = 0 ) = B^o 
p(x, y,z;t = 0) = po exp 


exp 


Lg 


y 




exp - 




( 1 ) 


where oJcs(x) shows the characteristic width at different height of the initial current sheet 

and Lg =120 Mm is the scale hight for a fully ionized hydrogen plasma: 




d ' x^ x-\- xq 

X Xq 


and Bxo, B^q, po, d, xq are normalized quantities: B^q = 0 . 2 , po = 1 . 0 , Bq 
d = 0.003 and xq = 20.0. 


( 2 ) 

= 1 - 0 . 


The X and y components of the magnetic field (5^, By) are obtained from the magnetic vector 
potential A as B = V x A. Note that the magnetic field strength slightly decreases via oJcs(x) 
with height ’x’ corresponding to the magnetic field in the solar corona balancing the gravity force 
(Eq.([^). The initial magnetic field state is displayed in the left panel of Fig[^ 


Compressible, resistive MHD equations (Eqs. 0 to ([^) are solved to describe the evolution of 
the plasma and magnetic fields (e.g., |Priest|1984| ): 


^ + V • (pu) = 0 

Ot 

dn 

p-^ +p(u-V)u = -Vp + JxB+pg 

Ot 


5B 

dt 


-V X E = V X (u X B - 77 J) 


dU ^ ^ 

^ + v.s = p«.g 


(3) 

(4) 

(5) 

( 6 ) 


where p is the plasma density, u plasma velocity, B magnetic field strength, E electric field strength, 
rj resistivity, g gravitational acceleration at the photospheric level and p plasma pressure. The 

Article number, page 6 of[^ 













X. Zhou et aL: Electron acceleration by cascading reconnection in the solar corona 


current density J, total energy density U and energy flux S are deflned as: 


VxB 


J = 

Po 

u = 



y - 1 2 

s = 



pu + 

2/io 


u B 


B + —J X B 

J^o j^o 


(7) 

( 8 ) 
(9) 


where yo = ^ is the adiabatic coefficient for adiabatic condition and po is the vacuum magnetic 
permeability. The anomalous resistivity j] in the Ohm’s law (Eq.([^) and in the energy flux S 
(Eq.(|^) is chosen ad hoc to describe the sub-grid-scale dissipation effects of microphysical (ki¬ 
netic) processes. It is switched on depending on the strength of the local current-carrier drift ve¬ 
locity vccD = \ J\/(ep) compared to the critical threshold velocity Vcr (e.g., Barta et al.|2Qll| ): 


T](r, t) = 


0 

^ vccpiX. t)-Vcr 

Vo 


IvcCdI < Vcr 
IvcCdI > Vcr 


( 10 ) 


Buchner & Elkina|2005[|2006| and |Karlicky & Barta|2008| have confirmed this behaviour by means 


of Vlasov and PIC-code numerical simulation and derived both critical velocity Vcr and the coeffi¬ 
cient C in Eq.(p^. 



Fig. 1. Evolution of the in-plane magnetic field companies of cascading magnetic reconnection in the CME- 
trailing current sheet obtained by high resolution 2.5D AMR MHD simulation. Panels from left to right show 
the initial state (t = 0 to), primary plasmoids (t = 80 to), secondary plasmoids (t = 200 to), third stage of 
plasmoids (t = 300 to), large scale magnetic islands mature state (t = 360 to) and last state (t = 520 to) where 
the erupted and disconnected magnetic field lines imply the appearance of a CME. 


The AMR MHD simulation (Eqs.([^ to ^) is carried out with normalized parameters: the 
normalizing length scale (half width of the current sheet at x = 0) is chosen to be Lq = 6.0 x 10^ 
m, the normalizing magnetic field is = 4.0 x 10“^ T and the normalizing number density is 
tzo = 1.25x 10^^ as well as = \e\ is taken as the normalizing charge. Other scaling parameters 
can be derived as: Vq = Bq/ = 7.80 x 10^m/s (where mo = Mp- proton mass and Vq is the 

asymptotic value of the Alfven velocity at y ^ oo, x = 0 and = 0), The time is normalized by the 
Alfven transit time to = Lo/Vo = 1.69 x 10“^ 5'. Eurthermore there are Eo = VoBo = 3.12x lO^V/m 
and qo = PoLoVo = 5.88 x lO^Q • m. The asymptotic plasma beta parameter isyS = 0.1 at (y ^ oo 
, X = 0). In the coarse resolution, the mesh sizes are Ax = Ay = 0.045 Lq. Hence the whole 
simulation domain extends over (0,288) x (-36, 36) in the x-y plane. 
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Fig[2 depicts the evolution of the in-plane magnetic fields. Total simulation is performed over 
520 to when a CME is ejected through the upper boundary of the box (last panel of Fig[^. We 
pick out an already fragmented current sheet at r = 360 as the background electromagnetic fields 
since at this time step not only there are the most information of the refined smaller-scale magnetic 
structures but also after that no more additional magnetic islands are generated. No anomalous 
resistivity is switched on before t = 420 ^o- 

In order to relate the electron acceleration to the resolution of the magnetic structures, we 
compare the acceleration in the coarsely and finely resolved magnetic fields. Fig [^compares the 
magnetic structures obtained by the coarse (upper panels) and higher (lower panels) resolutions at 
t = 360 ^ 0 - From left to right, increasing zoom-levels show the details of the magnetic structures. 
The right bottom panel depict the detail of smaller-scale magnetic structures obtained by the higher 
resolution. 



-5 0 5 -4 -2 0 2 4 -0.3-0.2-0.10.0 0.1 0.2 0.3 


Y/U Y/L„ Y/Lo 

Fig. 2. Magnetic field lines at time t - 360 to- Top line: coarsely resolved magnetic structures. Bottom line: 
higher resolution simulation. Left to right: increasing zooms. 


3. Methods Used 


3.1. Test Particle Calculations 


mV 2.71171 

If the gyroradius (rgy = —) and gyroperiod (oc 1 /ojgy = -) of the particle are much smaller than 

qB qB 

the length scale of transverse gradients (r^) and characteristic oscillation periods (oc 1/oJos) of the 


ambient electromagnetic fields (i.c., VgyIr_i_ 1 and ojgylojos » 1), a guiding center approximation 
is valid. The motion of a magnetized charged particle can be decomposed into a drift of its guiding 
center and a gyration around this center ( |Northrop|1963| ). 

The minimum magnetic field strength obtained by the AMR MHD simulations is 0.19 for 
10 MeV energized electrons, the corresponding gyroradius is 4.4 m only. The grid size even of the 
refined mesh (Ax = Ay = 0.0045 Lq = 2.7 km) is much larger. As well as, in normalized Eqs.(pT|) 
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to with the normalization values shown in Sectj^ a coefficient 
in Eqs.(p^ and ( p3] ). 


moVo 

qoBoLo 


and its reciprocal arise 


corresponds to the ratio of the particle gyroradius over the 


characteristic length Lq or the particle gyro-period 
mo Vo 


mp 

^ 0^0 


to the scaling time tp = LpjVp. If only 


is much smaller than unity, the guiding center approximation can be applied. In our study, 

qpBpLp 

it is only of the order of 10“^. Hence, here we use the guiding center approximation to trace each 


electron. 


Although only 0.01% and 0.45% electrons can be accelerated up to energies > 100 keV, for a 
high precision, a relativistic guiding center approximation is used: 


^ ^ (rvii) 

= vd+ b 

at y 

M (yv\\)^ m u 

Vd = Ve + - X (b • V)b] +- ^[b x (V(fcB))] 

q yk^B q yk^B 

(11) 

+ -^[b X (b • V)ve] + - ^[b X (ve • V)b] 
q yk^B q yk^B 

my 1 Em 

(12) 

‘^*^'11) =‘^ E - b ^[b-V(tB)] 

at m y 


+ (rv||)vE • [(b • V)b] + yvE • [(ve • V)b] 

(13) 

ic’^ + (yvii)^ + 2iaB 

(14) 

dt 

(15) 


here R, Vd, vy, y and b are the guiding center position vector, the perpendicular drift velocity, the 

c 

) and the magnetic field direction 


velocity along the magnetic field, the relativistic factor ( 

B 




unity vector b = —, respectively. In the expression for the drift velocity vd in Eq.JT^, the term ve 
B 

E X B 

corresponds to the local E x B drift velocity Ve = —^. Other terms are the magnetic curvature 
drift velocity and the magnetic gradient drift velocity as well as higher order drifts. The factor 

k = Jl 


Ve 


relates the electromagnetic field values to the reference frame moving with the 

If 

velocity Ve- Finally, yu = —- is the relativistic magnetic moment per mass unit where 

aB 

is the particle gyration velocity perpendicular to B. The electron energy is expressed using the 
relativistic y-factor as E = (y - l)mc^. The set of Eqs. O to ([T3 are solved utilizing a fourth- 
order Runge-Kutta scheme. The field values between the grid points are interpolated along the 
electron trajectories with 2D linear interpolation. 


4.752 X 10^ test electrons are initially uniformly distributed along the current sheet (0 < v < 108 
Lo, y = 0) at 2400 points with 22 different initial velocities from 0.0 to 21.0 Vth and 9 different initial 
pitch angles from 0 to tt. Here Vth is the electron thermal velocity for a typical coronal temperature 
of 10^ K\ 0.76Vo = 6 X 10^ km/s. Every electron is traced for up to lO^o (~ 0.769 s) or until it 
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leaves the simulation domain, whatever happens first. Not that this time is shorter than the time 
scale of essential magnetic field changes in the MHD simulations. 

3.2. Spectrum Distribution function of accelerated electrons 

To obtain the energetic electron distribution function, we use the fact that the solar corona is practi¬ 
cally collisionless. Hence according to Lioville’s theorem, the particle distribution function keeps 
constant along the particle trajectory: f(E,A, r, t) = f(EQ,A(), fq, ^o)- This allows to calculate the 
electron distribution function f(E, A,r,t) at the place where HXR are expected to be generated by 
Bremsstrahlung of energetic electrons. 


3.3. HXR Emission 


Knowing the local plasma number density and electron distribution function, the hard X-ray emis- 
sivity 1(6) integrated over all contributing electrons can be calculated in the frame work of the thin 
target model ( |Brown|1971| ) as: 

oo 

/(e) = ^ n(r)v(r)o-B(e, E)fiE, A, r, t) (16) 

E>e 


Here A, r and v(r) are the energy, pitch angle, position and velocity of the electron at the time t, 
n(r) is the local plasma number density, e is the radiated photon energy, f(E, A, r, t) is the electron 
distribution function at the position of interest place and time, while ctb is the cross section of the 
Bremsstrahlung process. For a simple approximation, we take the Bethe-Heitler formula for the 
Bremsstrahlung cross section ( |Bethe & Heitler|1934t Brown|197T ): 


(Ts(6, E)oc --In 
eE 



(17) 


Note that the Bethe-Heitler formula applies only to particle energies less than 100 keV. In this 
investigation, for both kinds of magnetic fields resolution, more than 99% of the electrons are 
accelerated to energies less than 100 keV, i.e. the Bethe-Heitler formula still can give a high 
accuracy here. 


4, Results 

Depending on the locations of the simulated electrons at 10 to, three groups of electrons can be 
identified: those trapped in the magnetic islands; those precipitating to the chromosphere and the 
ones being ejected into the interplanetary space. There is no electron escaping from the left and 
right sides of the simulation domain. We concentrate our analysis on the trapped (Sect |4.1| ) and 
precipitating (Sect |4.2| ) electrons. 
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4.1. Trapped Electrons 


More than 80% of simulated electrons are still trapped in the magnetic islands along the current 
sheet by 10 ^o- This highly dynamical and complex magnetic field structure provides a very effec¬ 
tive trapping mechanism of energetic electrons for the coronal HXR sources. 


4.1.1. Acceleration dependence on initial conditions and magnetic field resolution 

The acceleration of electrons in the convective (or induced) electric field E = -u x B is sensitive 
to the initial position, velocity, and pitch angle of injected electrons and the fine structure of the 
magnetic field. The upper panels of Fig[^ depict the dependence of the energy gain on the initial 
conditions and the magnetic field resolution. The lower panels show the corresponding projected 
results. 

In general, the electron acceleration is more efficient in magnetic fields with better resolved 
small scale structures for the larger magnetic curvatures and gradients accessible. The maximum 
final kinetic energy of trapped electrons is at most of the order of 100 in coarsely resolved 
magnetic fields, but it can be up to 470 keV if smaller-scale magnetic structures are taken into 
account, corresponding to a maximum energy gain of 53 keV and 420 keV for the coarsely and 
finely resolved fields respectively. 



Fig. 3. Dependence of the electron kinetic energy gain (AE) on the initial pitch angle (Aq), velocity (vq) and 
position (xo) for two differently resolved magnetic field structures. Each point represents one electron. The 
upper panels show the 3D result with the initial velocity color coded. The lower panels show the correspond¬ 
ing projected results and the averaged value (★ lines) and standard deviation of the energy gain (☆ lines) in 
the initial energy and pitch angle spaces. By along the current sheet center is depicted in the bottom-right 
panel. The red and blue colors are used to distinguish the results in the coarsely and finely resolved magnetic 
fields, respectively. Note that there are three different scales in the y-axis for AE < 0, 0 < AE < 3 keV and 
AE>3 keV. 
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From left to right the bottom panels of Fig [^depict the dependence of acceleration efficiency on 
the electron initial energy, pitch angle and position, respectively. The kinetic energy gain increases 
with the increase of the initial energy, which is consistent with the results of Guo et aL]|2010[ It 
is interesting to note that both the mean and the standard deviation of the energy gain are roughly 
proportional to the initial energy and the acceleration efficiency of the finely resolved case is about 
3 times higher than the coarse one. The dependence of the electron energy change on the initial 
pitch angle and position, however, is more or less chaotic due to the complex field structures. Dif¬ 
ferent from |Karlicky & Kosugi|2004| where the betatron process dominates, here the most energetic 
electron is not associated with an initial pitch angle of 90° any more. The acceleration symmetry 
with respective to the 90° pitch angle is also broken when the magnetic fields are better resolved. 
The lower-right panel of Fig also shows the magnetic field component By along the current sheet 
center y = 0, which can be used to identify the magnetic X- and 0-points with By = 0. The most 
efficient acceleration appears to be associated with electrons injected close to the X-points that 
contain larger magnetic gradients and smaller magnetic curvature radii. 


4.1.2. Energy gain 

The guiding center approach decomposes particle energy into components parallel and perpen¬ 
dicular to the magnetic field and the part associated with the guiding center drift in the direction 
perpendicular to the magnetic field. The maximum drift velocity Vu (Eq.(p^) is 1.40 Vth and 0.97 
Vth in the coarsely and finely resolved magnetic fields, which is negligible comparing with the other 
two components. 

Considering that the anomalous resistivity is not switch on (E • b = 0) and vd = Ve, ^ = 1, 
Eqs.([^, ([^ and ^ give: 

= - AiV||[b • VB] + (yviif ve • [(b • V)b] (18) 

dfiB dB 
dt ^ dt 

=/dV\\[h • VB] yUVE • (19) 

So the energy evolution of an electron in the guiding-center limit is given by: 
dEk fi?(rv||)^ + {yv^_f 

- CX- 

dt dt 

=n\^ ■ VB + (7 V||)^ve • [(b • V)b] (20) 


1 ^(yv||)^ 

2 dt 

1 djyvj^ 

2 dt 


Fig|^ exhibits the spacial distributions of the acceleration rates in Eqs.lp^ and ( p^ : (b • 
VB)I{2B) (the middle panel), (ve • VB)I{2B) (the left panel) and Ve • [(b • V)b] (the right panel) 
in the coarse calculation domain (corresponding to the top-middle panel of Fig|^, note that these 
distributions are quite similar between the coarsely and finely resolved magnetic fields. 
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Eqs.lp^ and ( p^ depict the parallel magnetic gradient (h VB) can change both the parallel 
and perpendicular energies of electrons, but they cancel out each other in the total energy evolution 
( |20| ). The overall energy change is dominated by the perpendicular magnetic gradient (jj\e • V5) 
and curvature (7V||)^Ve • [(b • V)b], which are proportional to the electron energy due to yu and 
(v||)^, respectively. The increase of the electron acceleration with the increasing initial energy is 
therefore expected (see the lower-left panel of Fig|^. (In addition, trapped electrons with larger 
initial velocities bounce and pass accelerators more frequently to gain more energies.) 


Ve_gradientB_360to 


Parallel_gradientB_360to Ve_curvatureB_360t|, 



Fig. 4. Spacial distributions of perpendicular (left panel), parallel (middle panel) gradient and perpendicular 
curvature ((right panel)) in Eqs.ip^ and ( p^ of the coarsely resolved magnetic fields. 


Also due to the combined actions between the magnetic gradient and curvature in Eq.(|^, the 
favourable initial pitch angles (shown in the lower-middle panel of Eig[^ are not 0, 180° or 90° 
which correspond to acceleration dominated only by magnetic curvatures ((7V||)^Ve • [(b • V)b]) or 
gradients (/iVE-V5), respectively. No favourable initial pitch angle with 0, 180° or 90° indicates that 
the perpendicular magnetic gradient and curvature acceleration efficiencies are comparable with 
each other in this complex magnetic field structure no matter what resolution of magnetic fields 
are used. These favourable initial pitch angles also change with the magnetic field resolution. The 
acceleration asymmetry around initial pitch 90° (the bottom-middle panel of Eigj^ is due to the 
non-symmetric acceleration factors around the current sheet center in (jjy^ • VB) and (7V||)^Ve • [(b • 
V)b] by the third dimension of electromagnetic fields in 2.5D symmetric current sheet geometry. 
The acceleration symmetry is weakly broken for the coarsely resolved case. 



Fig. 5. Parallel (AFy in top line) and perpendicular (AE^ in bottom) acceleration symmetry about the initial 
pitch angle 90°. Left panels: case of the coarsely resolved magnetic fields. Right panels: highly resolved ones. 
Every trapped electron (shown by one point) is color-coded by its initial velocity. The corresponding color- 
code is shown in the middle. 
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Fig|^ shows electron acceleration symmetry around the initial pitch angle 90° for the paral¬ 
lel and perpendicular energy gain components. With the details of the parallel and perpendicular 
acceleration, one can see that electron parallel and perpendicular acceleration are not exactly sym¬ 
metric around the initial pitch angle 90° in both the coarsely and finely resolved magnetic fields. 
The symmetry is better preserved in the coarsely resolved fields due to the smoothing effects (see 
the bottom-middle panel of Fig|^. 

The term yUV||(b • WB) in Eqs.(p^ and ( p^ , however, is not influenced by the third dimension of 
electromagnetic fields, hence it is symmetric around the current sheet center at the beginning. In 
other words, the non-symmetric acceleration in the parallel and perpendicular direction (see the top 
and bottom panels of Fig|^ respectively) are due to non-symmetric 0 uve-V 5) and ( 7 V||)^VE-[(b-V)b] 
around the current sheet center, respectively. For the coarsely resolved case, the acceleration is 
dominated by the perpendicular component, the reverse is true for the finely resolved case. Mean¬ 
while non-symmetric parallel acceleration (i.e. non-symmetric |v|||) can enhance the asymmetry in 
both the parallel and perpendicular acceleration (see Eqs.lp^ and (p^). 


E,,/Vs. Coarse E,,, Vs. E,^ Fine 



Fig. 6. Comparison between the final parallel (F||) and perpendicular (Ej_) kinetic energy of trapped electrons 
in the coarsely (left panel) and finely (right panel) resolved magnetic fields. Each electron is color-coded by its 
total kinetic energy change (AE). The black lines in each panel correspond to’E^ = Fy ’ and 'E - E^ +E\\ - 50 
keV\ 

Fig 1^ shows the distribution of the trapped electron acceleration in the (E\\, E’^) plane. The 
distribution of the strongly accelerated electrons is highly anisotropic. It is dominated by the per¬ 
pendicular and parallel kinetic energy components for the coarsely and finely resolved magnetic 
fields, respectively. While for the weakly accelerated electrons, they still roughly keep their initial 
isotropic distribution. For electrons initially moving only along magnetic field lines (with an initial 
pitch angle 0 or 180°), they do not have acceleration in the perpendicular direction. In this condi¬ 
tion, however, the parallel acceleration is a little stronger in the coarsely resolved magnetic fields 
(maximum finale energy 81 keV) than that in the fine case with maximum finale energy 79 keV. 

4.1.3. Characteristic trajectories 

To better understand the details of the electron acceleration processes, the first row of Fig|7]and 
second row of Figj^ depict the trajectory and energy evolution of the most energetic electrons in 
the coarsely (vq = 21 Vth, = 5/8;:, vq = 97.56 Lq) and finely (vq = 21 Vth, Aq = 31471, vq = 41.58 
Lq) resolved magnetic fields, respectively, and all panels in Fighave the same initial conditions 
as the corresponding panels in Fig|7] 
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Fig. 7. Trajectory and energy evolution for three characteristic trapped electrons in the coarsely resolved 
magnetic fields. Panels in the^r^t and second column show the electron trajectory in the vy and yz plane, they 
are color-coded by the total kinetic energy profile. Panels in the third column show the electron total kinetic 
{black line), parallel {red line) and perpendicular {blue line) energy evolution. 


For the coarsely resolved case, the most efficient acceleration is dominated by the increase of 
the perpendicular energy component for the positive perpendicular magnetic gradient //Ve • 
above 95 Lq (see the left panel of Fig|^. Also the slightly decreasing parallel energy is due to 
the negative perpendicular curvature Ve • [(b • V)b] there (see the right panel of Fig|^. With the 
same initial conditions, the corresponding electron in the finely resolved magnetic fields gain more 
energy (see the first row in Fig[^. 

For the finely resolved case, the most efficient acceleration happens to the parallel energy com¬ 
ponent with slightly increased perpendicular energy when the electron are trapped in the magnetic 
island (around x = 37 Lq) and accelerated again and again with its circulating motions by the posi¬ 
tive perpendicular magnetic curvature Ve • [(b • V)b] (located at the thin layer in the central current 
sheet above x = 40 Lq, corresponding to the step-like increased displacement in the z direction, 
see the final kinetic energy color-coding yz-trajectory projection in the second row of Fig[^ and 
gradient //ve • (located around x = 40 Lq). As expected, in the coarsely refined magnetic fields 
without the smaller-scale magnetic field structures, trajectory and acceleration of the corresponding 
electron are totally changed. This electron is mirror-trapped and cannot circulate in the magnetic 
island and move systematically in the z direction. 

In the right panel of Fig|^ one can see the largest perpendicular curvature acceleration region 
is located at around x = 85 Lq, however, the strongest energetic electron is not launched there. That 
is due to the cancellation between the perpendicular curvature acceleration and deceleration when 
electron circulates in the magnetic island at x ~ 90 Lq. Electrons in the third row of Fig|7]and Fig[^ 

Article number, page 15 of[^ 















X. Zhou et aL: Electron acceleration by cascading reconnection in the solar corona 


are chosen to reveal the acceleration characteristics around x = 90 Lq. These two electrons also 
have the same initial conditions: vq = 21 v^/^, Aq = 7/8;: and xq = 96.075 Lq. Their color-coded yz- 
trajectory projection in the third row of Fig|7]and Figj^prove the above discussion and their energy 
gains are mainly due to the parallel acceleration by the perpendicular magnetic curvature. Further¬ 
more these two electrons have the same trajectory and energy evolution, i.e., there is no influence 
coming from the magnetic field resolution, since no refined smaller-scale magnetic structures are 
found along their trajectories between x = 81 - 97 Lq. 



Fig. 8. Trajectory and energy evolution of electrons launched with the same initial conditions as the corre¬ 
sponding one in Fig|^but in finely resolved magnetic fields. 


Energy oscillation between parallel and perpendicular energies in each characteristic electron 
energy evolution profile (the last column of Fig|7]and Fig[^) is due to the parallel magnetic gradient 
V||(b • VB) in Eqs.(p^ and ( p^ when electron passes the positive and negative parallel magnetic 
gradient regions in turn or electron is mirrored with alternate parallel velocity in the parallel and 
anti-parallel direction. Each condition can be found in Fig|^and Fig While the magnitude of this 
oscillations is due to magnitude of the parallel magnetic gradients b • along electron trajectory 
(see the middle panel of Fig|^. 


4.1.4. Comparison with Observations 


More than 60% of trapped electrons are accelerated (AE > 0) and more than 50% of them have ki¬ 
netic energies larger than 10 keV. These energetic electrons can produce HXRs by Bremsstrahlung 
(note that the HXR range is 10 - 400 keV). In the framework of the thin target model ( |Brown| 
\91\) and using the Bethe-Heitler formula for the Bremsstrahlung cross section ( Bethe & Heitler| 
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see details in Sec |3.3| ), we derive the HXR spectrum of these energetic trapped 
electrons in order to compare our results with solar flare HXR observations. 

Since the initial electron distribution function in the solar atmosphere is not known, we consider 
three different initial distribution functions as: 


T9^|Brown|T97T 


/(£’o,Ao,ro,ro) ^ 


E’q (or Constant) 

Maxwell-Boltzmann 


( 21 ) 



- Coarse : 1 ; 11 

— ■ — - - ■ Coarse Weighted E"® : 4 ; 12 

— — — Coarse Weighted_Maxwell : *** ; 
- Fine : 1 ; 4 

-- Fine Weighted_E“® : 4 ; 5 

— — — Fine Weighted Mcixwell : *** ; * 



Energy (keV) 

- Coarse : 2 ; 12 

■ ■ ^ ■ ■ ■ Coarse Weighted E”® : 5 ; 13 

— — — Coarse Weighted Maxwell : *** 
- Fine : 2 ; 5 

■ ■ ■ ■ ■ Fine Weighted_E“® : 5 ; 6 

— — — Fine Weighted Maxwell *** ; 


Fig. 9. Electron {top) and HXR {bottom) spectra of energetic trapped electrons in the coarsely {red lines) and 
finely {blue lines) resolved magnetic fields with three different initial distribution functions - constant {solid 
lines), power-law -3 {dash-dot lines) and Maxwellian at 10^ K {dashed lines). The spectral indices are for 
the ranges marked with the black dashed lines embraced by two plus signs at two ends and their values are 
shown under each panel: the first one for electron and photon energies below 50 keV and the second one for 
the energies between 50 keV and 100 keV. 

The resulting electron and HXR spectra and spectral indices (below and above 50 keV) after 
acceleration {t = 10 ^o) are depicted in the top and bottom two panels of Fig|^ respectively. In 
approximation, the relationship between the electron {ye) and corresponding HXR {juxe) spectral 
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indices agree well with the relationship yuxR = + 1 in the thin target model. The influence from 

the ambient plasma number density Hy (Eq.(p^) is very small due to it small normalized range 
[0.3 - 1.8] in both differently resolved magnetic fields, while the HXR flux at 100 keV differs by 
more than 4 orders of magnitude for these two cases (see the bottom panel of Fig|^. Note that 
the HXR spectral indices, calculated from an initial Maxwell-Boltzmann distribution function for 
T = 10^ K, are too large to match any observed HXR spectrum. For the cases with a power-law 
distribution initially (e.g., [Karlicky & Barta|[2006| ), one may treat the electron acceleration as a 
diffusion in 2D energy space. Since the diffusion coefficient is approximately proportional to the 
energy, this explains the difference by one of the spectral indexes of the injected and accelerated 
electrons below the maximum injection energy of ~ 50 Above 50 the acceleration in the 
finely resolved magnetic field is much more efficient than that in the coarsely resolved magnetic 
field, we have a harder spectrum (~ 6) for the finely resolved case. The corresponding HXR spectral 
indexes are consistent with the observed values for small flares (whose HXR spectral indices can 
be as soft as > 7, see | Aschwanden|20'02| ) . 


EndPos Vs. E5„,| 0 Coarse 


EndPos Vs. E5„3| O Fine 




Fig. 10. Final locations of trapped electrons with final kinetic energies > 10 keV at t = lOto- They are color- 
coded by their final kinetic energies. Left and right panels correspond to the coarsely and finely resolved 
magnetic fields. Note that the electrons with final kinetic energy > lOSkeV are shown only by red asterisk 
points in the better resolved magnetic fields. 


Besides the HXR spectra, fine structures (bright spots) along the current sheets trailing CMEs 
or eruptive filaments were observed (e.g., by [Ciaravella et al.||2002[ [Ko et al.||2003[ [Savage et al. 
2010| ). These bright spots should come from energetic trapped electrons. Fig[T0| shows the final 
locations of these trapped electrons with final kinetic energies > \^keV which will brighten the 
magnetic island in the current sheet that may be associated with the observed hot spots. Further¬ 
more depending on the evolution of these magnetic islands, the bright spot located at v = 90 Lq 
moves upwards away from the sun while others fall back to the sun. This evolution agrees well 
with the observed upward ( [Ciaravella et al.|2Q'Q2{|Ko et~ar 2003}[Savage et al. 2010 1 and downward 
( [Ko et al. [2003[[Savage et al.|2010[ ) moving bright spots in CMF-trailing current sheets. 


4.2. Precipitating Electrons 

A second observable feature which can be derived from our calculations are the emission produced 
by the energetic precipitating electrons. They can precipitate to the solar chromosphere and be 
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related to the observed footpoint HXR signatures there. We first study the acceleration dependence 
of precipitating electrons on the initial conditions (velocity, pitch angle and position). The results 
are shown in FigfTT] 



* standard Deviation 


Fig. 11. Same as Figj^but for precipitating electrons. Here three different scales in the y-axis are also used 
for AE < 0, 0 < AE < 2 keV and ^E > 2 keV. 


4.2.1. Initial condition dependence 

Similar to electrons trapped in the current sheet, acceleration of precipitating electrons also strongly 
depends on their initial (velocity, pitch angle and position) conditions. The acceleration efficiency 
increases with the increase of the energy and the overall acceleration is more efficient in finely 
resolved magnetic field. However the acceleration is much less efficient than those trapped elec¬ 
trons. The maximum energy gain is only a few keV and about 10 for the coarsely and finely 
resolved magnetic fields respectively. The dependence of the energy gain on the initial pitch angle 
and position show that only electrons in a few channels can escape from the acceleration site and 
injected into the chromosphere. As expected, electrons moving along magnetic field line are more 
likely to escape than those with a pitch angle close to 90°. However, only a small portion (< 12%) 
of electrons can precipitate into the chromosphere. 

The bottom right panel of Fig [TT] depicts that a large portion of precipitating electrons start near 
to X-points. But no electron escapes from the X-points near v = 42 Lq and v = 97 Lq. Since 
the magnetic islands below this two X-points is not symmetric about their center - 0-points: their 
upper parts are smaller than the lower parts (see the whole By plots in the bottom right panel of 
Fig|^. Current sheet center launched electrons are easily reflected or trapped by this geography 
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Fig. 12. Same as Fig|^but for precipitating electrons. 


(see the characteristic trajectories of trapped electrons launched near x = 42 Lq in the middle-line 
of Figs |7] and [8]). 


4.2.2. Acceleration properties 

Although there are the same reasons for the asymmetric acceleration around the initial pitch angle 
90^ in the finely resolved magnetic fields between trapped and precipitating electrons, acceleration 
asymmetry of precipitating electrons is much weaker than that of trapped electrons (comparing 
Fig 1^ with Fig[T^. 

The total acceleration of precipitating electrons also is much weaker than that of trapped elec¬ 
trons, see the = 50 keV parts of Figand Fig[^ especially the coarse case for precipitating 
electrons in Fig[^ The final kinetic energy Ee of the most energetic precipitating electrons is a 
little more than 50 and 60 keV in the coarsely and finely resolved magnetic fields, respectively, 
i.e., all precipitating electrons have final kinetic energies < 100 keV. Different from trapped 
electrons in Fig|^ here most precipitating electrons still keep their initial energies shown as stripes 
parallel to = 50 

Acceleration difference between trapped and precipitating electrons is mainly contributed by 
the acceleration in the parallel direction. Precipitating electrons have stronger deceleration than 
acceleration in the parallel direction (see the top panels of Figand Fig[T^, i.e., the acceleration 
of precipitating electrons are mainly coming from the perpendicular direction independent on the 
magnetic field resolution (see also Fig[T^. For a stronger parallel acceleration, electron should stay 
longer around the current sheet center where has larger magnetic curvatures than other places, while 
precipitating electrons are ejected out of the current sheet before they can reach higher energies. 



Fig. 13. Same as Fig[^ but for precipitating electrons. 
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On the whole the acceleration differences of precipitating electrons between in the coarsely 
and finely resolved magnetic fields are not so many as that of trapped electrons, since precipitating 
electrons spend most of their time far away from the central current sheet and the better resolved 
magnetic structures are located only near the current sheet center. 


4.2.3. characteristic trajectory of precipitating electron 



Fig. 14. Trajectory projection and energy evolution of the most energetic precipitating electron with initial 
velocity vq = 21 Vth and pitch angle 157.5° in the coarsely {left panel) and finely {right panel) resolved 
magnetic fields. Trajectories are color-coded according to the local electron total kinetic energy. 


Top and bottom lines of Figshow the trajectory and energy evolution of the most energetic 
precipitating electron in the coarsely and finely resolved magnetic fields, respectively. Their vy 
and yz trajectory projections (along the magnetic field lines only) prove their strongly magnetized 
condition. Also the energy profiles in the last column of Fig[^ depict the different acceleration 
properties in the coarsely and finely resolved magnetic fields: their acceleration sites are still in the 
current sheet center by the perpendicular magnetic curvatures Ve • [(b • V)b] and gradients Ve • V5. 
After they leave there is no acceleration any more, there parallel magnetic gradient V||(b • VB) is 
stronger than the other two terms (see Fig|^. 

Also because of the single sign of the parallel magnetic gradients and direction of parallel 
velocity along precipitating electron trajectory, precipitating electrons do not have frequent energy 
oscillation as that of the characteristic trapped electrons in Fig |7] and Fig[^ 


4.2.4. Comparison with UV and EUV observations 

The low energies of the precipitating electrons in the convective electric fields can not cause HXR 
emissions but ribbons of UV and EUV brightening ( [Fletcher et al.|2Ql ij ). Fig[T^depicts the spatial 
distribution of the electrons precipitating to the chromosphere at the end of calculation (t=10 1^, top 
panels) and their evolution with time (panels in the last two lines). As the figure shows the ribbons 
exhibit a anti-symmetric geometry around the PIL. 

This two ribbons are related to the initial pitch angles of precipitating electrons: electron with 
an initial pitch angles > 90° (< 90°) precipitates into one (the other) branch. This kind of initial 
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Fig. 15. Chromosphere locations of precipitating electrons at t=10 to {top line) and their evolution {bottom 
two lines), color-coded by their final kinetic energies {E^ - Blue for < 54 keV and red for E^ > 54 
keV) separately in the coarsely and finely resolved magnetic islands. 


pitch angle dependence is attributed to the weak parallel accelerations by the perpendicular mag¬ 
netic curvature (7V||)^Ve • [(b • V)b] (Eq.(p^) which cannot accelerate electron into the direction 
anti-parallel to its initial velocity. Because of the non-symmetric acceleration around the initial 
pitch angle 90° of precipitating electrons in the finely resolved magnetic fields, more efficiently 
accelerated precipitating electrons are with initial pitch angles 0.875;: (> 90°). Hence more accel¬ 
erated {Ee > 54 keV) precipitating electrons are only located at one branch of the ribbon geometry 
with finely resolved smaller-scale magnetic structures (the top-right panel of Fig[T5]). Some of the 
observed asymmetry between two footpoints therefore may be attributed to the acceleration pro¬ 
cess. While the chromosphere energy distribution of precipitating energized electrons, accelerated 
by coarsely resolved magnetic structures, are more anti-symmetric with respect to the PIL (the 
top-left panel in Fig[T5]). 

Also with the chromospheric location evolution of the precipitating electrons (panels in the 
last two lines of Fig[T5]), one can find their locations along the chromospheric ribbons depend on 
their initial positions also: electrons started closer to the sun surface precipitate closer to the PIF 
, earlier in the chromosphere and have shorter displacements along z-axis (or PIL). At the same 
initial position, electrons with larger initial energies correspond to larger final kinetic energies and 
parallel velocities which lead electron to reach chromosphere earlier (see FigfT^. 

Fig[T^ shows that fluxes of electrons with higher energies evolve faster and reach peaks earlier 
than those of lower energy electrons. The time scale of the flux peak of precipitating electrons with 
final kinetic energies > 5QkeV indicates precipitating electrons are accelerated less than 1.0 < 0.1 

V. The refined magnetic field structures are mainly located above x = 25 Lq, hence at the beginning 
(before 1.5 ^o) there is no acceleration difference between the coarsely and finely resolved magnetic 
fields for precipitating electrons. 
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Fig. 16. Lightcurve of precipitating electrons for four energy ranges: Ee < 10 keV - dashed lines 10 < < 

25 keV - dotted lines 25 < Eg < 50 keV - solid lines Eg > 50 keV - dash-dot lines for the acceleration in the 
coarsely {red lines) and finely {blue lines) resolved magnetic fields. 

5. Conclusions and Discussion 

5.1. Conclusions 

In contrast to acceleration in direct current (DC) parallel electric fields which in MHD simulations 
depends on the choice of the resistivity in the Ohms law, we concentrate on the acceleration due 
to magnetic gradient and curvature drift effects in the cascading reconnection current sheet. We 
found that both trapped in the magnetic islands and precipitating electrons can be accelerated by 
the perpendicular magnetic gradients and curvatures. Trapped energetic electrons contribute to the 
formation of bright spots along the current sheet trailing CMEs or eruptive filaments) as well as the 
flare loop-top HXR radiation by their Bremsstrahlung. Precipitating electrons, on the other hand, 
cause ribbons of UV and EUV brightening in the solar chromosphere. 

Whether an electron becomes trapped or precipitating depends on the initial conditions (e.g., 
for precipitation, an electron should have a position around X-points, velocity > 2 Vth and pitch 
angle 9 ^ 90°). Trapped electrons are energized mainly in the magnetic islands in the coarse mag¬ 
netic fields, while in the better resolved magnetic fields, the strongest trapped electron acceleration 
takes place close to the X-points due to there finely resolved larger magnetic curvatures and gradi¬ 
ents in the smaller-scale magnetic fields. Both trapped and precipitating electrons are accelerated 
or decelerated in dependence on their initial positions and pitch angles. The electron final kinetic 
energy strength depends on the initial electron energy - larger initial energies cause stronger accel¬ 
erations. As well as every kind of electron can get more energization if the smaller-scale magnetic 
structures, obtained by higher resolution MHD simulations, are taken into account. Also because 
of these smaller-scale structures, energization of more accelerated trapped electrons are mainly in 
the parallel direction in the finer magnetic fields. Other (less accelerated trapped and precipitating) 
electrons mainly gain energies in the perpendicular direction. Due to the asymmetry in the mag¬ 
netic curvature drift acceleration term around the center of 2.5D current sheet, the larger magnetic 
curvatures in the better resolved magnetic structures cause stronger non-symmetric accelerations 
around initial pitch angle 90° of trapped and precipitating electrons. On the contrary, in the coarsely 
resolved magnetic fields both trapped and precipitating electron acceleration are close to symmetric 
around 90°. 
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With the better resolved small-scale magnetic structures, the maximum energy gain of trapped 
electrons can be up to 421 keV. This already suffices to explain the observed loop-top HXR ra¬ 
diations. Under the thin target model together with a simple Bethe-Heitler formula for the cross 
section of Bremsstrahlung and an initial distributions function cx (or constant) and oc E~^ , the 
HXR spectral indices of trapped electrons can be as hard as ~ 5 in the better resolved magnetic 
fields. This is already hard enough to explain the observed HXR spectra in medium solar fiares. 
For initial Maxwell-Boltzmann distributions for T = 10^ the HXR spectra provide, however, 
just a slight enhancement of the high energy tail. 

In the chromospheric ribbon-shape locations of precipitating electrons, electrons starting lower 
in the solar atmosphere precipitate closer to the PIT. The weak parallel acceleration of precipitating 
electron leads electrons with initial pitch angles < 90° precipitate to one side of the PIL, while ones 
with initial pitch angles > 90° go to the other side of the PIL. Generally, there is a anti-symmetrical 
geometry of precipitating electron locations in chromosphere around the PIL. While because of 
the stronger accelerations of electrons with initial pitch angles > 90°, more energetic electrons are 
located in one side of the PIL only with the better resolved smaller-scale magnetic structures. 


5.2. Discussion 


Solar flare observations imply that a large number of energetic electrons should precipitate into the 
solar chromosphere where they cause observable radiations. Our calculations have shown that only 
12% electrons can precipitate within 10 ^o- While the whole current sheet evolution is as long as 
520 to. Depending on the magnetic field evolutions (see panels of Fig[^, lower magnetic islands 
(x < 70 Lo) in Fig[T0|) will merge into one magnetic loop (see right panel of Fig[^ eventually. So 
in the end, the electrons previously trapped in the lower magnetic islands can later also precipitate 
to the chromosphere. Taking into account this merging effects, more than 63% electrons will 
finally reach the chromosphere. As well as when the space scale collapses to the kinetic one, the 
guiding center approximation will be not valid any more. Particle motion will become chaotic 
due to nonlinear resonances between particle bounce motion and gyration. With the transition to 
chaos, [Buechner & Zelenyi|1989| found that trapped nonadiabatic charged particles can escape due 
to chaotic pitch angle scattering effects. Furthermore in this study during 10 the background 
electromagnetic fields are constant, so the time effects on electron acceleration are neglected. With 
the evolutions of the electromagnetic fields, maybe some trapped electrons become precipitating 
ones. As a result, even more electrons will precipitate. 


Our study can explain the observed medium and small solar flare loop-top HXR spectra and 
EUV-ribbons just based on magnetic gradient and curvature effects in magnetic islands without 
ad hoc postulated "anomalous" resistivity. Precipitating electrons in our results, however, cannot 
explain the HXR spectral indices in the foot-points of solar fiares which can be as hard as 1.5 
in large solar fiares. Precipitating electrons also could reach the energies necessary to explain 
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the observed HXR spectra bases on magnetic gradient and curvature effects if more and smaller 
magnetic field structures are formed by cascading magnetic reconnection. 


Appendix A: Validation of the accuracy by using the conservation of the 
Second Adiabatic Invariant 


The conservation of the second adiabatic invariant ( [Northrop]|1963| ) of trapped electrons can be 
used to validate the accuracy of the numerical scheme solving Eqs. O to ([T§: 


J\\= m Vjjdl = constant (A.l) 


In Eq.( |A.l| ), the integral is taken along the particle guiding center trajectory between the mirror 
points ’a’ and ’b’. Eig |A.l| shows an example electron with conserved /||. 




. 

TTl 





> 


t/t„ 


Fig. A.l. Characteristic electron orbit indicating the conservation of the second adiabatic invariant. The 
upper left panel shows the XY-projection of the trajectory and the blue * points in this panel corresponding to 
the mirror points. The upper right panel shows /|| values along the trajectory and its red, green and deep pink 
* points show absolute changes of J\\ during half-period (deep pink * points) or one-period (red and green * 
points). Three kinds of velocity (perpendicular gyration velocity - blue line, parallel velocity - red line and 
drift velocity - dark line ) and energy (perpendicular gyration energy - blue line, parallel energy - red line and 
total kinetic energy - dark line) are separately shown in the bottom left and right panel. 


The upper left panel of Eig |A.l| depicts the XY-projection of the electron trajectory. Every 
blue asterisks in the upper right panel of Eig |A. 1 1 corresponds to a mirror point. Note that in the 
definition of the second adiabatic invariant /|| (Eq.( |A.l| )), dl > 0 when the particle is moving 
forward and dl < 0 when it is moving backward. Eor a convenience here we do not change the 
sign of ’dr when calculating the values of J\\ along the electron trajectory. One easily obtains 
not only the constant /|| over half-period (deep pink points) but also the vanishing J\\ over a 
whole bounce period (red and green points) as one can see in the upper right panel of Eig |A.l[ 
A comparison of parallel and perpendicular gyration velocities with the drift velocities along the 
electron trajectory (bottom left panel of Eig |A.l| ) illustrates the validation of the second adiabatic 
invariant. The bottom right panel of Eig |A.l| indicates that the total kinetic energy is exchanged 
between the parallel and the perpendicular directed motion. 
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